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Abstract 



The contribution of A hyperons to neutrino scattering rates is calculated in the random 
phase approximation in a model where the interaction is described by a Skyrme potential. 
Finite temperature and neutrino trapping are taken into account in view of applications to 
^Nq' ' the deleptonization stage of protoneutron star cooling. The hyperons can remove the problem 

of ferromagnetic instability common to (nearly) all Skyrme parametrizations of the nucleon- 
nucleon interaction. As a consequence, there is not any longer a pole at the transition in the 
^ , neutrino-baryon cross section. However there still remains an enhancement in this region. In 

the absence of ferromagnetism the mean free path in nph. matter is reduced compared to its 
, value in np matter as consequence of the presence of this additional degree of freedom. At high 

' density the results are very sensitive to the choice of the A-A interaction. 

o ■ 

^ . PACS numbers: 26.60. +c, 25.30.Pt, 21.60.Jz 

I keywords: neutrinos, proto-neutron stars, non-relativistic RPA, Skyrme potential, hyperons 

^ . 1 Introduction 
^3 . 

P5 , During and shortly after the collapse of a supernova, a copious amount of neutrinos is produced. 
The hope that these neutrinos could permit a successful explosion in numerical simulations through 
the Wilson (delayed) mechanism, together with the increase of available computational power, has 
triggered an important effort in designing a new generation of computer codes. The Boltzmann 
\ equation for neutrino transport can now be solved simultaneously with the equations for the hydro- 
dynamical evolution of the supernova explosion (see e.g. [1]). The differential cross section for the 
various neutrino production, absorption and scattering processes is a key input in the Boltzmann 
equation. These calculations also yield the spectral characteristics of the corresponding neutrino 
burst, which are now well within the reach of modern detectors. Although the current picture is 
that neutrinos alone are not enough to prevent shock stall, but 3-D fluid flow with convection is 
also necessary, the solution of the Boltzmann equation still seems to be an unavoidable requirement 
of a supernova explosion simulation. It is therefore necessary to have as good as possible a descrip- 
tion of the neutrino interaction with the matter and of the influence of medium effects, chemical 
composition and temperature on this data. 

The neutrino-nucleon cross section and the influence of nuclear correlations on this parameter 
in the central protoneutron star have been estimated by several authors both in relativistic and 
nonrelativistic models [2]-[13]. While it is generally found in relativistic models that medium 
effects tend to reduce the scattering rate, the findings from nonrelativistic models are less clear 
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cut. The nonrclativistic calculations on the other hand have been performed with a greater variety 
of assumptions as to the nature of nuclear correlations taken into consideration as the relativistic 
ones. 

A key issue is whether a ferromagnetic instability is excited at high density. In this case it 
would appear as a pole in the neutrino-nucleon cross section when the nuclear correlations are 
calculated in the random phase approximation [10, 11]. Such a feature is commonly observed when 
the nuclear interaction is described by a Skyrme force. It has been traced to the exchange (Fock) 
term. A ferromagnetic transition is actually also possible in relativistic models [18, 19] treated 
in the Hartree-Fock approximation. A relativistic calculation of neutrino-nucleon scattering at 
the RPA level and including exchange terms however has yet to be performed. Then again in 
nonrelativistic models the susceptibility has been found to increase and a ferromagnetic state to 
be energetically unfavourable in recent calculations where the nuclear interaction is extracted from 
Brueckner-Hartree-Fock calculations [12, 13]. 

At later times hyperons begin to appear and affect the luminosities as well as the equation 
of state [14]. One can hope to detect the tail of this neutrino burst (between 20 and 50 seconds 
after the beginning of the signal), which corresponds to the cooling of the still hot but partially 
deleptonized protoneutron star. In these later stages the diffusion equation is thought to be a good 
approximation, the neutrino dynamics then enter in the mean free path. After ~ 50 s, the neutron 
star has cooled to temperatures of the order of 1 MeV and achieved full deleptonization, i.e. the 
neutrinos are no longer trapped. It then enters in a long cooling phase. The mean free path is no 
longer the relevant parameter; instead, the available cooling codes need as an input the neutrino 
emissivity. 

The aim of this work is to assess the relevance of the hyperonic degree of freedom on the 
neutrino interaction with neutron star matter with a simple description of the interactions, before 
appealing to the heavier machinery of microscopical description through the Brueckner-Hartree- 
Fock calculations or taking into account the relativistic effects. We have chosen for this purpose 
the Skyrme interaction parametrized by Lanskoy el al. [15, 16, 17]. This description includes the 
A hyperons with parameters fitted with the available data on hypernuclei. It is also able to give a 
description of neutron stars in agreement with the results of more refined models [20] . 

Our model is described in section §2 with a presentation of the chosen Skyrme parametrizations 
in §2.1 - §2.3 and of the necessary formalism for the calculation of the cross section in §2.4. Section 
§2.5 discusses the influence of hyperons on a possible ferromagnetic instability and section §3 gives 
an estimate of the evolution of the hypcron fraction during the deleptonization of the proto-neutron 
star. We then present results for the vector and axial response functions and differential cross section 
in section 4.1 and for the mean free path in protoneutron star matter in section §4.2. Section §5 
summarizes the main results and gives a brief discussion of how the model could be improved. 

2 Formalism 

2.1 Skyrme forces 

Besides the usual parameterization of the nucleon-nucleon interaction 

VNNiri - ra) = to (1 + xoP„)5{ri - ra) + ^ti (1 + xiP^) [/c'^ S{ri - rs) + <5(ri - rs) 

+t2 (1 + X2Pa)k' 5{n -r2)k + (1 + x^Pa) p% S{n - r2), (1) 
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we use the following Lambda-nucleon and Lambda-Lambda potentials as suggested e.g. by Lanskoy 
et al [15, 16] 

VNAifN - ta) = -uo (1 + yoPa)S{rN - ^a) + ^ui fc'^ S{rN - ta) + S{rN - rA) 

+U2 k' 5{rN -rA)k + ^u^ (1 + ya^'a)^^ i^-^^-p^) '^(^^ - ^a) (2) 

VAA(ri - r2) = Ao(5(ri - ra) + ^Ai [fc'^ 5{ri - r^) + 5{ri - r^) k'^ 

+X2k' 5{ri - r2)k + XspAplf (3) 

(We haved dropped in these expressions the spin-orbit terms which will not be used in the remainder 
of this paper). 

While this model only includes the A hyperon, and would seem less complete than e.g. the model 
of Balberg and Gal [21] which includes all the hyperons, we fixed our choice on the former because 
these authors provide the two-particle potential. This permitted us to derive the Landau parameters 
for the spin 1 as well as for the spin channels. The Landau parameters in the spin 1 channel are 
namely necessary for the calculation of the axial structure function in the RPA approximation. By 
contrast, Balberg and Gal [21] give directly a parametrization of the (unpolarized) energy density 
functional. This parametrization can therefore only be used in the mean field approximation. While 
thermodynamical relationships can give access to the Landau parameters in the spin channel, we 
have no way of determining those of the spin 1 channel. Since the axial structure function is known 
to give the dominant contribution to the neutrino cross section, this precludes the use of the model 
of Balberg in the RPA approximation. We will nevertheless use the model of Balberg and Gal at 
the mean field level in order to estimate the error committed by neglecting other hyperons and in 
particular the S"~. 

Some a posteriori justification may be given for a model with the A hyperon only: (i) The 
parametrization of the NA and AA is relatively well grounded on Brueckner-Hartrcc-Fock calcula- 
tions and duly cross-checked with data on single and double A hypernuclei. By contrast, experi- 
mental data concerning the S hyperon is quite scarce, (ii) It is argued from recent analysis of H" 
atoms and the absence of a clear signal for the formation of bound S states in {K~,it) reactions 
[22, 23, 24], as well as from calculations in quark models [25, 26], that the S single particle potential 
is probably repulsive in nuclear matter. /3-equilibrium calculations performed on the basis of this 
assumption predict that the threshold for S hyperon formation in neutron stars is shifted to very 
large densities. If such be the case, they would concern only a minor fraction of the star or even 
not appear at all (see e.g. [21, 27]). (Hi) We are still concerned here with investigating qualitative 
effects, namely the influence of a new degree of freedom on the scattering rate and the effect of 
hyperon formation on preventing an eventual ferromagnetic transition. 



2.2 Composition of matter in chemical equilibrium 

In neutron star matter we impose that the conditions for /3-equilibrium are fulfilled. We write 
the equality of the chemical potentials fi = (pn + i^n) — [pp + 'nrip) = Pe — Pu, Pn + = 
PA+fnh- We must also impose electric charge conservation rig = Hp. At a given baryonic density the 
electron, proton and hyperon fractions are determined by the solution of these three equations. The 
chemical potentials are determined by deriving the Skyrme energy density functional with respect 
to the density of corresponding particle. Their explicit expression was given in a companion paper 

[20]. The electrons are relativistic and their chemical potential is given by p^ = kp^ + m'^ . 

In protoneutron star matter with trapped neutrinos we have p,y = (67r^ray)^^/^\ while in colder 
neutrino-free neutron star matter p^ = 0. The muons were not considered for simplicity, as their 
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presence is not expected to lead to qualitative changes in the results. The threshold density for 
A hyperon formation is determined by the condition n\{thr) = ^AikpA = 0) = fin + "m-n — t^k- 
The hyperons form around 2-3 rigat in neutrino-free matter, they appear at higher density when 
neutrinos are trapped in the star. 



2.3 Skyrme parametrizations of the hyperonic sector 

In [20] we have examined many combinations of parameters for the Skyrme forces among those 
available in the literature in the light of their suitability for the description of neutron star physics. 
We settled on a recent parametrization by Reinhard and Flocard [28] in the nucleon sector SkI3 
(or the very similar SkI5) together with the set YBZ6 for the nucleon-A interaction among those 
recommended by Lanskoy et al. [15, 17] as offering one of the best description of hypernuclear data. 
For the A-A interaction we then considered two options, either the SLL2 force of Lanskoy [16] or 
no interaction. Lanskoy fitted the binding excess to ASaa = —4.8 MeV, i.e. more atttractive than 
the presently accepted value of —1 MeV. The truth should hopefully lie between these two choices. 



model 
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ti 


t2 


ts 


xo 


Xl 


X2 


X3 


SLylO 


1/6 


2506.77 


430.98 


-304.95 


13826.41 


1.0398 


-0.6745 


-1.0 


1.6833 


Skis 


1/4 


-1762.88 


561.608 


-227.09 


8106.2 


0.3083 


-1.1722 


-1.0907 


1.2926 


sv 


1 


-1248.3 


970.6 


107.2 


0. 


-0.17 


0. 


0. 


1. 



Table I: Skyrme parameters for the NN interaction [28, 29, 30] 
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100.4 


79.60 


2000. 


-0.107 


0. 



model 


Ao Al 


SLL2 
no LL 


-437.7 240. 




Table II: Skyrme parameters for the NA and AA interactions [15, 16, 17] 

All expected properties for the neutron star structure were then recovered. Very importantly, 
the ferromagnetic transition common to (nearly) all Skyrme models is avoided when hyperons are 
present (see §2.5). 

Among the main results obtained with the SkI3-|-YBZ6-|-SLL2 parametrization, let us mention: 

- The equation of state is softened by the presence of the hyperon degrees of freeedom. As a 
consequence, the maximum mass (2.26 Mq) for a npe star is reduced to 1.64 Mq for a star with 
hyperons. 

- A hot star with trapped neutrinos was observed to be mctastable in the mass range [1.64 - 
1.88] Mq, as noticed in several earlier studies, see e.g. [31]. A protoneutron star in this mass range 
will eventually collapse to a black hole as it undergoes deleptonization. 

- The equation of state remains causal in the whole star when hyperons arc present. 

- The A hyperons appear at a threshold density 2.08 rigat in cold neutrino-frec matter. In 
protoneutron star matter with a lepton fraction Yl = 0.4, the formation of hyperons is delayed 
until 2.85 ngat- 

- The asymmetry energy = {l/2)d{E/A)/df3'^\i3^o increases rapidly with density with the 
SkI3 parametrization, allowing for high proton fractions. 

- There exists a non negligible hyperon fraction in the core of the protoneutron star before 
deleptonization is fully completed (see section §3). The tail of the supernova neutrino burst which 
corresponds to this stage could therefore be affected by the presence of hyperons. 

Another suitable set was the Skyrme Lyon SLylO [29] together with the LYI parametrization 
of the N-A force also favored by Lanskoy. The asymmetry energy increases with density but much 
slower than with the SkI3 force; as a consequence, the proton fraction remains rather low. With 
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hyperons the maximum mass of the neutron star was a lower but still acceptable 1.425 Mq. The 
set SLylO was preferred over the more widely used SLy4 and SLy7 sets because it was easier to 
keep it clear of the ferromagnetic pole. In other respects the properties of SLylO are very similar 
to that of the Sly4 and Sly7 sets. 

Finally we singled out the older SV force [30]. This set is rather atypical as it has = 0, 
its effective mass at saturation is also lower (m* = 0.38) than the commonly accepted values. 
Although its adequacy for reproducing the properties of nuclei and nuclear matter is somewhat 
less satisfactory than was the case for SkI3 or SLylO, it has the unique feature of not presenting a 
ferromagnetic transition in nuclear matter nor in npe matter in /3-equilibrium which makes it very 
valuable for comparison purposes. Its general characteristics (stiffness and maximum neutron star 
mass, proton content) are quite similar to those of the SkI3 set. 

We reproduce below the relevant entries of the tables published in [20] 



Model 






^ferro 


SLylO+LY-I+SLL2 


2.719 


0.041 


6.183 


SkI3+YBZ6+SLL2 


2.076 


0.144 


9.194 


SV+YBZ6+SLL2 


1.951 


0.136 


9.093 



Table III: Thresholds for hyperon formation and for ferromagnetism 
in npAe matter in (3 equilibrium 
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nc(1.4 Mq) 


R{1A Mq) 


'^max 






SLylO, npe matter 


7.308 


3.60 


11.05 


7.69 


1.99 


9.52 


SLylO+LY-I+SLL2 


16.032 


9.38 


9.0 


12.38 


1.425 


8.11 


SkI3, npe matter 


6.343 


2.27 


13.21 


6.12 


2.263 


11.16 


SkI3+YBZ6, no LL 


causal 


2.36 


13.20 


4.66 


1.655 


12.52 


SkI3+YBZ6+SLL2 


13.077 


2.39 


13.19 


6.79 


1.642 


11.02 


SV, npe matter 


4.983 


2.10 


13.46 


5.44 


2.426 


11.54 


SV+YBZ6, no LL 


causal 


2.19 


13.40 


4.44 


1.662 


12.75 


SV+YBZ6+SLL2 


13.082 


2.22 


13.41 


6.5 


1.639 


11.24 



Table IV: Conditions for causality and neutron star properties 



2.4 Neutrino-beiryon scattering and absorption rates 

In the non relativistic approximation the differential scattering cross section for the process ^'(pi) + 
B{p2) —>■ i^iPs) + B{p4) is taken to be (see e.g. [2, 5, 10, 11]): 

= S^^^^'t^ ~ [(1 + cos0)5W(a;,fc) + (3 - cose)5«(c., fc)] (4) 

In this expression, lo and k are the transferred energy and momentum pj'—pg = (uj , k) , = E,^ — uj 
is the neutrino energy after the collision and 6 is the scattering angle pi-Ps = \pi \ {psl cos 9. The 
structure functions .S^*^) in spin channel S are defined by the expectation value of the density and 
spin density fluctuations in baryonic matter. In the mean field approximation they both reduce to 

= ^/(P2)(27r)^(5^(pi +P2-P3-P4) (1 - /(P4)) (5) 

They are related to the imaginary part of the reducible polarizations by: 

"^^^Ho^^fc) = -^YZr^^p^ IraU('\u,k) (6) 
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This definition differs by a factor l/vr from that of [10, 11, 32] but agrees with that of [4, 5] so that 
the elastic hmit be 50(0;, k) 2TmB5{ui) 

In the non relativistic limit, the spin and spin 1 responses correspond to the vector and axial 
couplings respectively. In the random phase approximation the vector and axial structure functions 
obey uncoupled Dyson equations: 



^vfA = no + ^vfA ^(5=0/1) no. 



where 



(27r) 



(7) 

(8) 



qo + Ep- Ep+q + ie 
and V(5^o/i) the interaction potential in the particle- hole channel. 

In our case, the Dyson equation has a 3x3 matrix structure. The RPA polarization is obtained 
by inverting (7) with 

no 
Hn = nH 





no 
no 



(9) 



We have for the real part of the polarizations 



with 



27r2A; 



pdp In 



- 2kpf - {2m*ujf 
(A;2 + 2kpY - (2m*a;)2 



/(^,) = |l + exp[(^-Ai)/r] + l 



and fn = iii 



(10) 



The expressions for the effective masses and chemical potential were given in the Appendix of [20] . 
At r = 0, TZe reduces to the Linhardt function in agreement with e.g. [32] 



7^e J^Q{uo,k) 



T-*0 



27r2 



(11) 



l + x 



1-x 



x± 



k , wm.* 



2pFi kpFi 



The imaginary part is always analytical 



Im U'Q{u,k) 



2iTk 



In < 



1 + exp 


/ P- 
\2m* 




1 + exp 


( pI 

\2m* 





with p± = PFiX± 



(12) 



We will extract the potential in the Landau Fermi Liquid approximation from the Skyrme 
parameterization. Our potential has no imaginary part in this approximation. We will be working 
in the simplest monopolar (/ = 0) Landau approximation; our potential has therefore no angular 
or momentum dependence. For a discussion of the full RPA with momentum dependence we refer 
the reader to [10, 11, 33, 34]. 

The potential in the vector isovcctor channel for neutrino scattering through the neutral current 
process is given in our model in terms of the Landau parameters fij 



V, 



(S=0) 



fnn fnp fnA 
fpn fpp fpA 
fhn fAp /aA 



(13) 
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For the axial channel, one should replace all fij by the gij in order to obtain the relevant potential 
^(5=1). The Landau parameters can be extracted from the Skyrme potential energy by applying 
a double functional differentiation with respect to occupation numbers [32, 33]. Their explicit 
expression when hyperons are present was given in [20] . 

After performing the inversion, we obtain for example in the axial {S = 1) channel 



n 



RPA 
Ann 



n 



RPA 

Anp 



n 



RPA 
A nA 



1 


Det[/- 




=1) 


no] 




1 






Det[I- 


Vis= 


=1) 


no] 




1 






Det[I- 




=1) 


no] 



(1 - gp^n^ii - saahI) - ^^An^no 

(1 - 9AAn.A)gpn + QpAgnA^A 
(1 - gppUl)gnA + gnp9pA^n 



^^p-^^n 



(14) 
(15) 
(16) 



with similar expressions for the vector channel {S = 0). These equations reduce to Eqs. (26-32) of 
Reddy et al. [7] when no hyperons are present. 

Let us note at this stage that if the determinant appearing in the denominator is strongly 
reduced or vanishes, this will give rise to a corresponding peak in the structure function and in the 
neutrino cross section. This is indeed what happens in the axial channel for previous calculations 
[10, 11] in pure neutron matter and npe matter in f3 equilibrium (see also Fig. 16 of [7]): The 
condition 



Da = Det[/ 



^(5=1) no 







(17) 



is fulfilled at some critical density, signalling the transition to a ferromagnetic state. This is 
discussed further in next section. 

Finally the structure functions appearing in equation (4) are obtained from 



with 
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<^A 



V/A '-V/A '-V/A 



) 



/ CPP CP" CP'i 

/ "^V/A "^V/A '^V/A 
'^V/A '^V/A 



V/A 



oAn 
•^V/A 



pA 

Oft A 
'^V/A 

oAA 
•^V/A- 



-(1 -4sin2 Ow) = 0.04 
1,^ 1.26 



1,^ 1.26 



f^/A 
^V/A 
V '^V/A , 



(18) 



= -^(F + L>/3) = -0.365 



(19) 



2.5 Poles in the dispersion relations 

Previous studies on the neutrino mean free path in neutron matter [10] or npe~ matter in (5 
equilibrium [11] found that a pole appears in the calculation of the axial structure function above 
a certain critical density. This feature is typical of Skyrme parameterization of nuclear interactions 
and is related to a transition to a ferromagnetic state. 

In [20] we defined the magnetic susceptibilities Xij where i,j G {n,;j. A}. The inverse suscepti- 
bilities are proportional to the second derivatives of the energy density functional with respect to 
the polarizations: 



Xij 



2p_ 
PiPj 



A, 



(20) 
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A, 



2 dsidsj 
Pi] - Pil 



(21) 
(22) 



Pi] + Pil 

On the other hand it can be shown that the Ajj are related to the Landau parameters through 

A _ 2p 1 



A,; 



2p 



r(l + Q) 



m*kFi 



(23) 



A criterion for the appearance of the ferromagnetic phase is that the determinant of the inverse 
susceptibihty matrix vanishes, or equivalently, in terms of the Landau parameters: 



Det 




G7 



(1 + Gin -0 





(1+oy 



(24) 



In the hmit where the temperature T and the energy transfer w go to zero, the real part of the 
Linhardt functions tend to the limiting value 



7^e^i)(u;,/^) 



{T->0, w->0, w/fc=CBt) 



(25) 



If we further assume that the imaginary part of the polarizations vanish (which, even though it 
turns out to be an excellent approximation, is not strictly the case) we would then find that the 
condition that a pole appears in the dressed axial polarization Det[/ — V^^^i^IIo] [cf. Eq. (17)] is 
identical to the criterion for the appearance of a ferromagnetic phase defined in this section. In 
fact the finite imaginary parts provide for some Landau damping of the pole. The pole will also be 
somewhat shifted and smoothed by finite temperature, but this hardly changes our conclusions in 
practice for the temperatures relevant for protoneutron stars (up to 30 MeV). 

In [20] we plotted the criterion (24) as a function of density. It was shown that when the 
threshold for hyperon formation was lower than the critical density for the ferromagnetic transition 
in npe matter, the hyperons were able to stabilize the system so that the pole could be avoided. 
An sample of these results is displayed in Fig. 1. 




Fig. 1 - Criterion for the transition to 
a ferromagnetic state. 
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A similar study can be performed for the vector response function. The denominator Dy = 
Det[/ — V(5=o)no] = reduces at T = 0, ^ to the criterion for instabihty under density 
fluctuations equivalent to Eq. (24) by replacing the Gij by the corresponding F^j. At low density, 
no hyperons are present, and the criterion reduces to that for spinodal instability in asymmetric 
nuclear matter (see e.g. [35]) {1 + Fq"'){1 + Fq^) — {Fq^)'^ = 0. The pole disappears at temperatures 
~ 15 MeV, but the enhancement persists at higher temperatures. 



3 Evolution of the hyperonic content in a deleptonizing proto neu- 
tron star 

As the protoneutron star cools and the trapped neutrinos arc released, the beta equilibrium becomes 
more favorable to the production of hyperons. If a sizable number of hyperons appear before 
the deleptonization is completed, the mean free path of the neutrinos will be affected and could 
influence the characteristics of the tail of the neutrino burst detected in supernova explosions. The 
determination of the evolution of the hyperonic content would require to solve self consistently 
the equations for thermodynamical conditions, neutrino mean free path, diffusion of neutrinos and 
cooling along the lines of reference [14] . 

In this section we take the results from Fig. 17 of Pons et al. [14] as representative of the 
evolution the temperature and leptonic content. With this input wc calculate the structure of the 
neutron star at each time step. While this procedure is not consistent, it still permits to estimate 
whether hyperons will be relevant at this stage. 

We chose the evolution curves correponding to model GM3 with hyperons (curved labeled 
GM3npH in [14]). The temperature T and total lepton number = Yg + Yu are fixed to values 
extracted from this figure (see table below). The central density was also taken according to the 
results of [14] but was rescaled for each parameter set in order to match with a 1.6 Mq neutron star 
at the end of deleptonization. It was checked that the equation of state and chemical composition 
at the end of deleptonization obtained by setting Yl = 0.09 coincides with that obtained when 
setting Fjy = 0. 
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Fig. 2 - Evolution of the hyperonic con- 
tent as a function of time 



After neutrinos have left the neutron star, it continues to cool from T = 18.3 MeV to T=0. We 
perfomed four more calculations at T=0, 5, 10 and 15 MeV to follow the hyperonic content in this 
later stage. The result is displayed in Figures 2 and 3 for the parameter sets which allow a stable 
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Fig. 3 - Particle fractions as a function of the 
neutron star radius at three stages of the delep- 
tonization. 
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neutron star with a mass 1.6 Mq, namely SkI3 or SV in the nucleon sector, YBZ6 for the hyperon- 
nucleon interaction and either no Lambda-Lambda interaction or the SLL2 parametrization of 
Lanskoy et al. 

In Fig. 2 the hyperonic fraction at the center of the star is displayed as a function of time. We 
can see that a sizable fraction of hyperons is present as early as 10 seconds after the collapse. The 
result are very similar whether we use the SkI3 or SV force. On the other hand we again note a 
strong dependence on the characteristics of the A-A interaction. 

Fig. 3 shows the density profile of the neutron star at three stages of the dclcptonization (t=20s, 
t=50s and t oo) for the parameter set SkI3-|-YBZ6-|-SLL2. As the star cools, its radius shrinks 
and the hyperons gather at its center. 

4 Neutrino scattering rates 

4.1 Structure functions 
a) Axial response 

Let us first discuss the behaviour of the contributions to the axial structure functions from 
neutrons {S^"'^), protons (S^^^) and A hyperons (S^^^). 
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As already obtained by other authors [10, 11], the axial structure function for neutrons in npe 
matter displays a strong enhancement at = 0, going over to a pole at the critical density ^if^j./o'^^- 
We have seen in §2.5, the criterion for ferromagnetic instability is equal to the condition that the 
denominator of the dressed polarizations vanish at a; = 0. This is a mechanical instability, meaning 
that an infinitesimal energy fluctuation uj = £ can trigger the phase transition. This should not be 
confused with the phenomenon of spin zero sound analogous to the zero sound mode observed in 
the vector reponse function, where a peak appears at a finite value co = lozs- In the latter case the 
response function displays a resonance corresponding to the excitation of a collective motion of the 
system at the corresponding velocity vzs- 

Fig. 4 compares the axial structure functions obtained in the mean field and random phase 
approximations. The system is npeA matter in P equilibrium at a density = 4 Usax, temperature 
T = 10 MeV and we fixed the transferred momentum at A; = 30 MeV. The interactions are described 
by the choice of parameters SkI3 + YBZ6 + SLL2. While the RPA correction is only moderate in 
the S^^^ and S^^^^ contributions, we can clearly see the enhancement due to the vicinity of the 

ferromagnetic criterion in S^"'^ . On the other hand, no spin zero sound peak is observed at a finite 
value of LO. 

Fig. 5 compares the S^"^ as a function of density for three different Skyrme parametriza- 
tions: SLyl0+LYI+SLL2, SkI3+YBZ6+SLL2 and SV+YBZ6+SLL2. At saturation the three 
parametrizations give very similar results. At three time saturation density on the other hand 
the quantity -Dyi(a; = 0) = det(l + Gij) is approaching the zero axis for the SkI3 and SLylO models 
(see Fig. 1) and the 5^"'' response is strongly enhanced. The SV+YBZ6+SLL2 parametrization 
which is free of ferromagnetism accordingly does not give rise to any enhancement at a; = 0. 

S^^^ follows the usual trend of having its strength shifted to higher oj when density is increased. 
This is also observed for S^"^ in the SV+YBZ6+SLL2 parametrization when the effect is not 
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Fig. 5 — Axial structure functions in 
npeh. matter in beta equilibrium in the 
RPA approximation - Comparison of 
three Skyrme parametrizations 



washed out by the enhancement at w = 0. At T = 5)^ ' \s only non zero above the threshold for 
hyperon production and then qualitatively behaves as S^^^\ At high density its contribution is of 
the same order of magnitude as that of the neutrons. 

Fig. 6 illustrates the temperature dependence of the axial structure functions. On the left panel, 
S^"^ is represented at = 4 ngat and fc=30 MeV. As the temperature is increased, the response 
function is enhanced and processes in which the neutrino gain energy in the collision {lo < 0) 
gradually open. On the right panel, the hyperon structure function is displayed at ub = 2 Usat, 
i.e. shortly below the threshold density n^^. = 2.08 ngat- While the structure function vanishes as 
expected at zero temperature, its contribution increases as temperature rises and rapidly reaches 
the same order of magnitude as scattering on nuclcons: Even though there are only traces of 
thermally created hyperons in this regime, they contribute actively to the scattering as they are 
not yet impeded by the Pauli blocking. 

b) Vector response 

In Fig. 7 wc show the vector structure function for the neutron. At low density we find when 
plotting Dy{(jj, k) = Det[/ — V(5=o)no] that this function comes near to the zero axis and eventually 
crosses it. This can be explained as the onset of the spinodal instability, when homogeneous matter 
is not the stablest state anymore, as is the case at the inner edge of the crust of the neutron star 
(see e.g. [36]). As the density is increased, this instability disappears but instead dips form in Dy 
for a finite value of the energy transfer. When the real part of Dy vanishes a zero sound mode 
develops. This occurs at approximately 3 psat depending on the equation of state. At moderate 
densities, the imaginary part is finite and provides for Landau damping of the mode. This is seen 
on the vector structure function as a peak with some width. At higher density the zero sound mode 
is unscreened and appears as a delta function outside the vector strength distribution. 

We observe that two such dips in Dy develop in npe matter in j3 equilibrium whereas a third dip 
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appears in npeK matter. This translates in respectively two or three peaks in the vector response 
function. The conditions of the right panel of Fig. 7 were chosen so that these features are clearly 
visible. It is tempting to ascribe each dip to one species of baryons present in the system. As 
a matter of fact, the position of the dips is approximately given by the characteristic frequencies 
associated to each particle tOi = k{k + 2pFi)/{2m*). The deviation from the exact value is due 
to the fact that we have mixed states: the isospin is not a good quantum number, neither is the 
strangeness. For the density range relevant to neutron stars we only see the zero sound associated 
to the neutrons. The position (i.e. velocity) of the zero sound mode(s) depends on the stiffness 
of the equation of state. On Fig. 7-b we see that, the equation of state with hypcrons being 
softer than in matter with only neutrons and protons, the zero sound mode propagates at smaller 
velocity. The trend is general, and the position of the zero sound peak for various choices of the 
NN, NA and AA forces follows their classification according to stiffness realized in [20]: SV-I-YBZ6 
> SkI3+YBZ6 > SLylO+LYI. 

Fig. 7-a was drawn at zero temperature. When some temperature is applied to the system 
(see Fig. 8), the delta function (corresponding to an undamped zero sound mode at T=0) merges 
with the remaining part of the structure function. As was already observed in the case of the axial 
response, processes in which the neutrino gain energy in the collision {u < 0) open as T increases. 

c) Differential cross section 

The differential cross section Eq. (4) is displayed on Figure 9 at = 2.5 rigat, *-e. shortly above 
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the threshold for hyperon formation. We compare the results for Skyrme models SkI3-|-YBZ6-|-SLL2 
and SV+YBZ6+SLL2. 

Tiny secondary peaks are barely visible at high energy transfer u; = it 60 MeV. This is what 
is left of the zero sound contribution in the vector structure functions (see Figures 7 and 8) after 
applying the Pauli blocking factor l — f(El). Moreover, in the chosen example Ei, = k = 30 MeV the 
condition | cos 0| < 1 restrict the allowed energy ransfer lo to the range [—k — 2 — k] = [—30 — 30] 
MeV, so that the zero sound peak does not contribute alltogether. 
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Fig. 9 - Differential cross section in npeA. matter in beta equilibrium in the RPA approximation 



4.2 Mean free path 

The mean free path is obtained by integration of the differential cross section. This parameter is 
relevant in the late stages of supernova explosion and first stages of proto neutron star cooling, 
when the neutrinos are still dynamically trapped in the hot and dense matter. 

^ dujdO. 



A J \ dwdVl 

Gl f^^ , E' f^E' 



47r2 



K / dcu-f / qdq [1 - f{El)] [(1 - cosd)Sv{cJ,q) + (3 + cos 9) SA{oJ,q)] 

r J-oo -c/i/ J\uj\ 

(26) 
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In the elastic limit where Sy ~ 



SI - 



27rnj(5(a;) for i,j G {n,p, A} and without Pauli 



blocking, this reduces to the textbook formula j = J2i 
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i=n,p,A TT K'^Vi ^'^Ai) ' 
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We studied the total cross section and the mean free path for both the neutrino-free and trapped 
neutrino case. According to protoneutron star cooling calculations [14], the characteristic energy 
of the neutrinos is ~ 3T in the former case whereas it is E^, ~ ji^, in the latter case. 
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Fig. 10 - Mean free path with and without hyperons in the mean field approximation (a) In 
neutrino free neutron star matter at T=5 MeV (h) In protoneutron star matter with trapped 
neutrinos at T=20 MeV, Yl=^A 



Fig. 10 shows the mean free path in the mean field approximation. On the left panel we 
have represented the neutrino-free case Yy = ^ for a moderate temperature T = 5 MeV while the 
right panel displays the case of trapped neutrinos 11 = 0.4 with a higher temperature T = 20 
MeV. It can be seen that, at the mean field approximation, the total cross section is enhanced by 
the apparition of hyperons and the mean free path correspondingly reduced. The same behavior 
was observed by Reddy et al. [6] and explained by the combined effect of the modification of the 
chemical potentials of the nucleons and the opening of new scattering channels when new degree 
of freedoms are available. The results of the SkI3-|-YBZ6-|-SLL2 and SV-I-YBZ6-I-SLL2 are similar, 
with the mean free path being larger for the SV parametrization corresponding to lower effective 
masses: At the mean field level, a oc Tmllo oc ml (sec Eq. (12)) and A rsj l/a. 

The mean free path is observed to increase with higher densities. This behavior was also found 
in nonrelativistic models by Reddy et al [7], whereas relativistic models yield a decreasing mean 
free path with increasing density. The discrepancy was explained by Reddy et al [7] as the result 
of two competing mechanisms, i. e. the decreasing of the effective masses (see previous alinca) and 
the increase of the chemical potentials: in the nonrelativistic Skyrme model the decreasing m* is 
the dominant effect. 

The mean free path in the random phase approximation is displayed in Fig. 11 in neutrino- free 
matter on the top panels and with trapped neutrinos on the bottom panels for parametrizations 
SkI3+YBZ6+SLL2 (left) and SV+YBZ6+SLL2 (right). 

We first note an enhancement of the cross section at subnuclear density for both parametriza- 
tions. This is related to the peak present in the vector response function which was noticed in the 
previous section, and signals the mechanical instability to non homogeneous "pasta" phases. Our 
model is not valid any more in this density range; a calculation would need to be performed along 
the lines of reference e.g. [37]. 

The dramatic enhancement of the cross section in the RPA approximation due to the onset 
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of ferromagnetism which was reported in several previous works [7, 10, 11] is also observed with 
parameter set SkI3 in npe matter. It is indicated on the figure by the line nfcrro- The mean free 
path in npe matter goes to zero at this point. When hyperons are present in the system, the 
ferromagnetic transition can be avoided when the threshold density for hyperon formation is lower 
than the critical density for ferromagnetism in npe matter [20]. For the set SkI3+YBZ6+SLL2 
we have n^j.esh = 2. 08 rigat while in npe matter nferro = 3.08 ngat- As a consequence, the pole is 
avoided and the mean free path goes on gently increasing. 

The parameter set SV which was selected as the only Skyrme model in which the ferromagnetic 
transition does not appear, permits us to study the influence of the RPA corrections in the absence 
of the pole. We see on the right panels of Fig. 11 that the mean free path behaves in this case in 
a way similar to the mean field result at densities np > ngat- 

Fig. 12 shows the ratio of the total cross sections URPA/frMF = Xmf / lambda^^pA in the RPA 
and mean field approximation. A residual enhancement in the vicinity of the ferromagnetic in- 
stability remaining from the np sector prior to the hyperon formation threshold is visible in the 
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SkI3+YBZ6+SLL2 parametrization. In the SV+YBZ6+SLL2 paramctrization the cross section is 
somewhat enhanced with respect to the mean field result: even in the absence of ferromagnetism, 
this is still at variance with the relativistic prediction [7, 8, 9]. 
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Fig. 12 - Enhancement factor for cross sections calculated in the random phase approximation 
with respect to mean field approximation (a) In neutrino free neutron star matter at r=5 MeV 
(h) In protoneutron star matter with trapped neutrinos at T=20 MeV, Yi,=0.4 

We finally explored the impact of the neglect of the other hyperons S, S by performing a mean 
field calculation in the model of Balberg and Gal [21] and another nonrelativistic model based on 

the Seyler-Blanchard effective interaction by Banik and Bandhyopadhyay [38]. Both models include 
all species of hyperons. The calculation was nevertheless performed here in a reduced version of 
these models with only the S~ and A hyperons taken into account and no muons. The S~ hyperons 
appear at 1.87 ngat and the A at 2.32 ngat in the model of Balberg and Gal with the parameter set 
corresponding to the choice 7 = 4/3, while the S~ are formed at 1.52 rigat and the A at 2.82 rigat 
in the model of Banik and Bandhyopadhyay. The mean free path obtained with such models is 
shown in fig. 13 in neutrino free matter at T = 10 MeV for a neutrino of energy Ey = 30 MeV . 
The reduction of the mean free path at the formation threshold of a new hyperon is also observed 
in these models. The effect is less visible for the second hyperon than for the first. The behavior 
of the mean free path obtained from the model of Banik and Bandhyopadhyay is qualitatively 
similar to that resulting from the Skyrme parametrizations considered in the present work. The 
decrease of the mean free path with increasing density obtained from the model of Balberg and 
Gal is easily understood from the fact that whereas both the Skyrme and the Seyler-Blanchard 
parametrizations have decreasing effective masses, Balberg and Gal keep theirs fixed at the free 
value (see the discussion below Fig. 10). 

5 Conclusion 

We have performed for the first time a calculation of the neutrino-baryon scattering rates including 
the hyperons at the RPA level in a nonrelativistic model. The nuclear interactions were described 
by Skyrme forces calibrated to reproduce the available data on nuclei and hypernuclei and tested 
for their suitability to describe the properties of neutron stars. 

At the mean field level the hyperons open new scattering channels and tend to decrease the 
mean free path as compared to its value in npe matter. In the random phase approximation the 
cross section is very sensitive to the degree of stability of the matter with respect to density or 
spin density fluctuations. The hyperons help delaying or removing the onset of a ferromagnetic 
instability usually present in Skyrme models. The drastic reduction of the mean free path at the 
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transition is thus avoided. A residual enhancement of the cross section nevertheless remains in 
most cases. 

Even in the complete absence of a ferromagnetic instability as provided by the parametrization 
SV+YBZ6+SLL2, we still obtained results which are at variance with the prediction of relativistic 
models: (i) The mean free path increases with density. This may be ascribed to a rapid decrease 
of the effective masses of the baryons in these models. This feature was observed previously in npe 
matter [7]. (ii) The RPA cross section was still slightly enhanced with respect to the mean field 
value. This indicates a less repulsive spin interaction as density increases. 

Our model is still very schematic. Several items would need to be improved in future work, 
among which: 

(i) Particle-hole interaction 
The Skyrme interaction was used for convenience in order to keep the model as simple as possible. 
While it was carefully selected in order to avoid known problems at high density such as the 
behavior of asymmetry energy, causality and onset of ferromagnetism, it is still being driven at the 
limit of its reliability. The behavior of the interaction in the spin-spin channel still remains unclear. 
An other issue is the description of the A-A force which plays an important role at high density 
where the hyperon fraction is large, whereas the description of the hyperon-hyperon interaction in 
phenomenological formalisms is still very poor. 

Both points could be substantially improved while still keeping the same level of simplicity for 
the Dyson equations and neutrino-baryon scattering cross sections, if we could obtain a parametriza- 
tion of the polarized energy density functional and effective masses from microscopical variational 
or Brueckner-Hartree-Fock calculations. Landau parameters extracted from BHF calculations have 
indeed been used very recently to calculate the neutrino scattering mean free path in neutron 
matter with the Nijmegen NSC97e potential [12] and for npe matter in f3 equilibrium with the 
Argonne AVig potential plus 3-body forces [13]. The spin-spin interaction was actually found to be 
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increasingly attractive by both groups, leading to a strongly reduced cross section in the random 
phase approximation. 

Brueckner-Hartree-Fock calculations also offer a coherent framework to calculate the hyperon- 
hyperon interaction in medium starting from bare baryon-baryon potentials and then cross checking 
the results with the data on double hypernuclei. The parametrization of Brueckner results for 
matter with an hyperonic component is currently under way and will be the subject of another 
paper [39]. 

(ii) Full RPA 

We limited ourselves to the monopolar I = Landau approximation. As previous studies in pure 
neutron matter or symmetric nuclear matter have shown [11, 33, 34] the full RPA approximation 
can represent an appreciable correction to the simple I = result. While the extension of the full 
RPA to npAe matter should not present any conceptual difficulties, the additional amount of work 
required is not yet worthwhile, considering the other uncertainties of the model. 
(Hi) Relativistic effects 

The main advantage of the nonrelativistic model is that the Dyson equation has a simpler 
structure, especially if we restrict the calculation to the monopolar approximation. This allows to 
explore a greater variety of many-body effects and chemical composition of the matter. At the high 
densities considered in protoneutron stars, a relativistic model however appears to be a desirable 
requirement. 

(iv) Charged current processes 

The present work explored the role of hyperons in neutrino processes on the example of the 
scattering through the neutral current. The same methods as used here can be applied to the 
charged current process. In protoneutron stars, the charged current process Vg + n ^ p + e~ 
is actually dominant. In the subsequent long cooling phase, the triangle rule for conservation of 
momentum constrained to the immediate vicinity of the Fermi surface is less restrictive in the 
direct URCA process involving hyperons such as p + e~ A + i^g, which makes this an interesting 
mechanism for an efficient cooling of the star [40] . 
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